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Abstract 

It is shown, in the framework of the Thomas-Fermi model at finite temperature, 
that a cooling non-degenerate gas of massive neutrinos will, at a certain temperature, 
become unstable and undergo a first-order phase transition in which quasi-degenerate 
supermassive neutrino stars are formed through gravitational collapse. For neutrinos 
in the mass range of 10 to 25 keV/c 2 , these compact dark objects could mimic the role 
of supermassive black holes that are reported to exist at the centres of galaxies and 
quasi-stellar objects. 

1 Introduction 

A gas of massive fermions, interacting only gravitationally, has interesting thermal properties 
that may have important consequences in astrophysics and cosmology. The canonical and 
grand-canonical ensembles for such a system have been shown to have a non-trivial thermo- 
dynamical limit Jl], |2|. Under certain conditions, these systems will become unstable and 
undergo a phase transition that is accompanied by gravitational collapse ||. It is interest- 
ing to note that this phase transition occurs only in the case of the attractive gravitational 
interaction of neutral particles obeying Fermi-Dirac statistics: it neither happens in the 
case of the repulsive Coulomb interaction of charged fermions [|J, nor does it in the case of 
the attractive gravitational interaction, when the particles obey Bose-Einstein or Maxwell- 
Boltzmann statistics. To be specific, we henceforth assume that this neutral fermion is the 
heaviest neutrino, which is presumably the z/ T , although this is not essential for most of the 
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subsequent discussion, as this neutral fermion could also be replaced with a light gravitino 
in supergravity models with low supersymmetry breaking scales ||37|| . We will, moreover, 
assume that the weak interaction can be neglected. 

The ground state of a gravitationally condensed neutrino cloud, with mass below the 
Oppenheimer-Volkoff (OV) limit where it is close to being a black hole, is a cold neutrino 
star H 0, ||, in which the degeneracy pressure balances the gravitational attraction of 
the neutrinos. Degenerate stars of neutrinos in the mass range between m u = 10 and 25 
keV/c 2 are particularly interesting [0, ||, jcTTR , as they could explain, without resorting to 



the black hole hypothesis, at least some of the features observed around the supermassive 
compact dark objects, which are reported to exist at the centres of a number of galaxies 
!, £3 0,0, including our own 0, [|| |7|, g|, M> H > and quasi-stellar objects 

(qso)ht 



18] , [19], ^0], ^9|. In fact, the upper bound for the mass of a neutrino star is given by 

where m Pl = (hc/G) 1 / 2 = 1.22205 -10 19 GeV/c 2 



pi 



m. 
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the OV limit M ov = 0.54195m 
is the Planck mass and g u the spin degeneracy factor of the neutrinos and antineutrinos 
(i.e. g u = 2 for Majorana or g v = 4 for Dirac neutrinos and antineutrinos). The radius of 
such a compact dark object would be i?ov = 4.4466-Rqv where Rq V = 2GM y/c 2 is the 
Schwarzschild radius of the mass M G v 
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There is thus little difference between a neutrino 
star at the OV limit and a supermassive black hole of the same mass, a few Schwarzschild 
radii away from the object, as the last stable orbit around a black hole has a radius of 
3 Schwarzschild radii anyway. For non-relativistic neutrino stars that are well below the 



OV limit, mass and radius scale as MR = 91.869/1 G 



3 m„ 8 
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For instance, the 



purported supermassive black hole Sgr A* at the centre of our galaxy could well be a neutrino 
star with 2.6 x 10 6 solar masses and a few tens of light-days radius |2S|, [51], [5jJ. 



The existence of a quasi-stable neutrino in this mass range is ruled out neither by par- 
ticle and nuclear physics experiments nor by direct astrophysical observations ||. On the 

v P oscillations [E2l, and the 



contrary, if the LSND experiment is actually observing — > 
quadratic see-saw mechanism involving the up, charm and top quarks, is the correct explana- 
tion for the smallness of the neutrino masses |23], Q , we should expect the v T mass between 
6 keV/c 2 and 32 keV/c 2 . In such a scenario, the v e and masses would be about 10~ 5 
and leV, respectively, and the solar and atmospheric neutrino anomalies would have to be 
interpreted as vacuum oscillations into two different sterile neutrinos, and z/jj , with masses 
slightly different from those of v e and z/^, consistent with the measured solar and atmospheric 
neutrino mass squared differences. Large neutrino-antineutrino asymmetries could suppress 
potentially disastrous active-sterile neutrino oscillations prior to nucleosynthesis |36 |. The 
mixing angles for active-sterile neutrino oscillations involving the same flavor would be max- 
imal, while the mixing matrix elements corresponding to flavor oscillations would be very 
small and exhibit a hierarchical pattern similar to that of the CKM matrix of the quark sec- 
tor. The smallness of the masses of the sterile neutrinos , and u~! could be guaranteed 
by an additional see-saw mechanism in the sterile sector p8| . 

It is well known that a neutrino mass of between 10 keV/c 2 and 25 keV/c 2 which is in 



the so-called cosmo logically forbidden region given by 93 h eV/c 



< 



m 



< 



where 0.4 ,5 h ^ 1 is the Hubble parameter, is unacceptable, as it would lead to an early 



< 



4 GeV/c 2 
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neutrino-matter dominated phase some time after nucleosynthesis and prior to recombination 
| 2j| . In such a universe, the microwave background temperature would be attained much 
too early to accommodate the oldest stars in globular clusters, nuclear cosmochronometry, 
and the Hubble expansion age, if the Standard Model of Cosmology is correct. However, 
the early universe might have evolved quite differently, in the presence of such a heavy 
neutrino. In particular, it is conceivable that primordial neutrino stars have been formed in 
local condensation processes during a gravitational phase transition that must have occurred 
some time after the nonrelativistic heavy neutrinos started to dominate this universe, a few 
weeks after the big bang. Aside from reheating the gaseous phase of the heavy neutrinos, the 
latent heat produced during the phase transition might have contributed partly to reheating 
the radiation as well. Moreover, the bulk part of the heavy neutrinos (and antineutrinos) 
will have annihilated efficiently into light neutrinos via the Z° in the dense interior of these 
supermassive neutrino stars || [?J ]8[]. In this context, it is interesting to note that, even 
within the (homogeneous) Standard Model of Cosmology, the reason why neutrino masses 
larger than 4 GeV/ c 2 are again cosmologically allowed, is precisely that the annihilation of the 
heavy neutrinos into light neutrinos reduces the matter content of the universe. Since both 
these processes, reheating and annihilation, will increase the age of the universe, or postpone 



the time when the universe reaches the current microwave background temperature [S, 21 



it does not seem excluded that a quasistable massive neutrino in the mass range between 10 
and 25 keV is compatible with cosmological observations f| |§. 



2 The Thomas- Fermi model 

The purpose of this paper is to study the formation of such a neutrino star as a consequence 
of a thermodynamical instability during a first-order gravitational phase transition. For 
simplicity, we assume that the neutrino star will be sufficiently below the OV limit, so that 
we can treat this process non-relativistically p5] . The general- relativistic extension of the 
Thomas- Fermi model has been discussed in [j30| , |33| . The effects of general relativity become 



important only if the total rest-mass of the system is close to the OV limit [51], [Jj| . The 
gravitational potential V(r) satisfies the Poisson equation 

AV = 47rGmln v , (1) 

where the number density of the r neutrinos (including antineutrinos) of mass m v can be 
expressed in terms of the Fermi-Dirac distribution at a finite temperature T as 



n v [r) 



9 » ,{2m v kTflH, (»- V W). ( 2 ) 



4vr 2 r v M kT 



Here I n (r]) is the Fermi function 



3 



and \i the chemical potential. It is convenient to introduce the normalized reduced potential 



' (/x-n (4) 



m u GM & 

M & being the solar mass, and the dimensionless variable x = r/Ro with the scale factor 

2/3 



Rn 



3tt/t 



±\/2mig v G*/ 2 M% 2 

{17.2 keV\ 8/3 2/ , 
2.1377 lyr ( - - j g' 2 ' 3 . (5) 



m,,c 2 



Using equations (2) and (4), equation (1) takes the simple form 



x dx 2 2 2 \ x 

where we have introduced the normalized inverse temperature f3 = T /T, with T = m u GM Q /kRo. 
In equation (6) we recover, at zero temperature, the well-known Lane-Emden differential 
equation || [8j 

d 2 v v 3/2 

5? = -j;- (7) 



The solution of the differential equation (6) requires boundary conditions. We assume 
here that the neutrino gas is enclosed in a spherical cavity of radius R corresponding to 
X\ = R/Rq. We further require the total neutrino mass to be M u , and we allow for the 
possibility of a pointlike mass Mb at the origin, which could be, e.g., a compact seed of 
baryonic matter. v(x) is then related to its derivative at x = x\ by 

1 / Mb + Mu\ , a , 
v (a?i) = — \v{x 1 ) — , (8) 



which, in turn, is related to the chemical potential by /i = kT v'(xi). v(x) at x = is given 
by the point mass at the origin, i.e. M B /M & = v(0). 

Similarly to the case of the Lane-Emden equation, it is easy to show that equation (6) 
has a scaling property: if v(x) is a solution of equation (6) at a temperature T and a cavity 
radius R, then v(x) = A 3 v(Ax), with A > 0, is also a solution at the temperature T = A A T 
and the cavity radius R = R/A. 

It is important to note that only those solutions that minimize the free energy are physical. 
The free-energy functional is defined as [0], 

F[n u ] = n[n v \N v - W[n v \ - kTg u J |^||x 
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In 



1 + exp 



p 2 V[n v ] n[n v ] 



1ra v kT kT 



+ 



kT 



(9) 



where 



and 



V[n v \ = -Gml J d 3 r 



\t 7*1 ' 



W[n u ] = —-Gmi I d 6 rd 6 r 



(10) 



The chemical potential in equation (9) varies with density, so that the number of neutrinos 
N v = M u /m v is kept fixed. 

All the relevant thermodynamical quantities, such as number density, pressure, free en- 
ergy, energy, and entropy, can be expressed in terms of v/x, i.e. 



M 3 



n,Ax 



r>% (^) 



M Q T 



x 



(12) 
(13) 



\liN v + \kT,RlJ d*xn v {x) v{x) x V{Q) 



Rl J d s xP„{x), 



(14) 



E = -fiN v - -kT R 3 J d 3 xn u (x 



x 



J d 3 xe kin (x), (15) 



S = ^(E-F). (16) 



3 Numerical results 



We now turn to the numerical study of a system of self-gravitating massive neutrinos with an 
arbitrarily chosen total mass M = 1OM , varying the cavity radius R. Owing to the scaling 
properties, the system may be rescaled to any physically interesting mass. For definiteness, 
the v T mass is chosen as m v = 17.2 keV/c 2 which is about the central value of the mass 
region between 10 and 25 keV/c 2 0, that is interesting for our scenario. In Figure 1, our 
results for a gas of neutrinos in a cavity of radius R = 100i?o are presented. We find three 
distinct solutions in the temperature interval T = (0.049 4- 0.311)Tq; of these only two are 
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Figure 1: Density distribution normalized to unity for condensate-like (solid line) and gas- 
like (dashed line) solutions at T — T t . Free energy, energy, and entropy per particle as a 
function of temperature. Temperature, energy, and free energy are in units of T . 

physical solutions, namely, those for which the free energy assumes a minimum. The density 
distributions corresponding to such two solutions are shown in the first plot in Figure 1. 

The solution that can be continuously extended to any temperature above the men- 
tioned interval is referred to as "gas", whereas the solution that continues to exist at low 
temperatures and eventually becomes a degenerate Fermi gas at T = is referred to as 
"condensate". In Figure 1, we also plot various extensive thermodynamical quantities (per 
neutrino) as functions of the neutrino temperature. The phase transition takes place at a 
temperature T t , where the free energy of the gas and of the condensate become equal. The 
transition temperature T t = 0.19442T is indicated by the dotted line in the free-energy plot. 
The top dashed curve in the same plot corresponds to the unphysical solution. At T = T t the 
energy and the entropy exhibit a discontinuity, and thus there will be a substantial release 
of latent heat during the phase transition. An important and currently still open question 
is, how and to which type of matter or radiation this latent heat, which can be interpreted 
as the binding energy of the neutrino stars, will be transferred. 

We acknowledge useful discussions with Dr. D. Tsiklauri and Dr. G.I. Kyrchev. 
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